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Abstract. Given the difficulty in directly determining prominence physical param- 
eters from observations, prominence seismology stands as an alternative method to 
probe the nature of these structures. We show recent examples of the application of 
magnetohydrodynamic (MHD) seismology techniques to infer physical parameters in 
C"| . prominence plasmas. They are based on the application of inversion techniques using 

Q_j observed periods, damping times, and plasma flow speeds of prominence thread os- 

cillations. The contribution of Hinode to the subject has been of central importance. 
We show an example based on data obtained with Hinode's Solar Optical Telescope. 
Observations show an active region limb prominence, composed by a myriad of thin 
horizontal threads that flow following a path parallel to the photosphere and display 
synchronous vertical oscillations. The coexistence of waves and flows can be firmly 
established. By making use of an interpretation based on transverse MHD kink oscil- 
lations, a seismological analysis of this event is performed. It is shown that the combi- 
nation of high quality Hinode observations and proper theoretical models allows flows 
and waves to become two useful characteristics for our understanding of the nature of 
solar prominences. 



1. Introduction 

Solar prominences are one the most intriguing manifestations of solar activity. These 
structures consist of large clouds of plasma, two orders of magnitude cooler and denser 
than the surrounding corona, suspended against gravity by forces thought to be of mag- 
netic origin. The physical properties of prominence plasmas are akin to those of the 
chromospheric plasma, hence some as yet not well determined mechanisms must pro- 
vide the required thermal isolation and hydrodynamic support during lifetimes that last 
from days to weeks. The nature of s olar prominence s is closely linked to their sub- 
resolu tion structuring. Early studies by de Jager ( 1959) and Kuperus & Tandberg-Hanssen 



( 19671) already pointed out that prominences are composed by many fine threads. This 



has been confir med by more recent high-resolution observations obtained by, e.g., 
Lin et al. I (l2005h . The fine thre ads are made of cool absorbing material, believed to 



outline magnetic flux tubes (Ma rtin et al.ll2008h . Their average width is about 074, their 
length is in between 5" and 40", and they have lifetimes of up to 20 minutes. 

The measurement of physical properties of prominence plasmas by direct obser- 
vation is challenging. An alternative to gain knowledge about the physical conditions 
in prominences is the combined use of observed and theoretical wave properties. The 
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Figure 1 . Systematic of MHD seismology. 



presence o f waves and oscillations in so lar prominences is known since a long time (see 
reviews by ICU lver 120091 : lBalleste3l201Ch . Individual threads or groups of them oscillate 



with periods that range between 3 and 20 minutes (lYi & Engvold 1991 :lYi et all 1991, 
LinetalJl2007l : lOkamoto et aDl2007l : iTerradas et al.ll2008l : iNing et al.H2009l : iLin et al. 



2009) . A recurrently observed property of prominence oscillations is their rapid tem- 



poral dampin g, with perturbations decaying in time-scales of only a few oscillatory 
periods (see lTerradas et alJl2002l : lLir]|2004l : lNing et al.ll2009l. for recent examples). An- 
other r elevant characteristic of waves in s olar prominences is the ubiquitous presence of 
flows dZirker et alJI 19981 : ILin et al.ll2003ll2005h . Flow speeds in the range 5-25 km s _1 
are usually reported in quiescent filame nt threads, while in active region threads flow 
speeds of up to 50 km s _1 are detected ( Ok amoto et al.ll2007h . 

Transverse thread oscillations are commonly interpreted in terms of standing or 
propagating MHD kink waves. The MHD wave interpretation of thre ad oscillations 
has a llowed the development of theoretical models (see reviews by iBallested [2005. 



2006). Theoretical d amping mechanisms have also been developed (see Ballester 2010; 



Arr egui & Ballestenl2010L for recent reviews). As information on geometrical proper 



ties, damping mechanisms, and flows are incorporated in theoretical models for MHD 
wave studies, more refined results are obtained that allow a better implementation of 
prominence seismology. 

Solar atmospheric seismology aims to determine physical parameters in magnetic 
and plasma structures that are difficult to measure by direct means. It is a remote di- 
agnostics method that combines observations of oscillations and waves in magnetic 
structures, together with theoretical results from the analysis of osc illat ory properties 
of mod els of those structures. It was first suggested bvlUchidal <jl970h and Robert s et al. 
(Il984l) . in the coronal context, and by Tandberg-Hanssen ( 1995) in the prominence con- 
text. The systematic of MHD seismology is displayed in Figured] Observations of solar 
coronal magnetic and plasma structures provide us with information from which theo- 
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retical models about their equilibrium can be constructed. On the other hand, observa- 
tions also provide us with measurements of certain properties, such as periods, damping 
times, or flow speeds. By analyzing the wave properties of given theoretical models 
they can be compared to the observed wave properties. If we find a perfect agreement 
between observed and theoretical wave properties we can aim to derive some unknown 
physical parameters of interest. Meanwhile, we can test and improve our models or find 
constraints. This article discusses three tools for the application of MHD prominence 
seismology, with an emphasis on the Hinode contribution to the area. 



2. Seismology Using Observed Periods of Thread Oscillations 

A recent application of the prominence seismology technique, using the perio d of ob- 
served filament thread transverse oscillations can be found in Lin et al. (2009). These 



authors find observational evidence about swaying motions of individual filament threads 
from high resolution observations obtained with the Swedish 1-m Solar Telescope in La 
Palma. The presence of waves propagating along individual threads was already evident 
in, e.g. jLin et al.l (120071) . However, the fact that line-of-sight oscillations are observed 



in prominences beyond the limb, as well as in filaments against the disk, suggests that 
the planes of the oscillation may ac quire various o rientations relative to the local solar 
reference system. For this reason, iLin et all d2009h combine simultaneous recordings 



of motions in the line of sight and in the plane of the sky, which leads to information 
about the orientation of the oscillatory plane in each case. Periodic oscillatory signals 
are obtained in a number of threads, that are next fitted to sine curves, from which the 
period and the ampl itude of the wave s are derived. The presence of different cuts along 
the structures allow ILin et al.l d2009h to obtain the phase difference between the fitted 



curves, which can be used to measure the phase velocities of the waves. The overall 
periods and mean velocity amplitudes that are obtained correspond to short period, P ~ 
3.6 minutes, and small amplitude ~ 2 km s" 1 oscillations. The information obtained 
from these Ha filtergrams in the plane of the sky is combined with Ha Dopplergrams, 
which allow to detect oscillations in the line-of-sight direction. By combining the ob- 
served oscillations in the two orthogonal directions the full vectors are derived, which 
sho w that the oscillat ory planes are close to the vertical. 

On et all (l2009h interpret the observed swaying thread oscillations as MHD kink 
waves supported by the thread body. By assuming the classic one-dimensional, straight, 
flux tube model a comparison between the observed wave properties and the theoretical 
prediction can be made in order to obtain the physical parameters of interest, namely 
the Alfven speed and the magnetic field strength in the studied threads. To this end 
the observed phase speed is directly associated to the kink speed, which in the limit of 
high density contrast, typical of filament plasmas, is simply reduced to c& V^v^,-, 
with VAi the prominence Alfven speed. This allows to obtain the thread Alfven speed 
through V Aj - Vphase/ V2. The obtained values for a set of 10 threads can be found in 
Table 2 in Lin et all i 2009h . Once the Alfven speed in each thread is determined, the 



magnetic field strength can be computed, if a given internal density is assumed. For a 
typical value of p,- = 5 x 10 -11 kg m~ 3 magnetic field strengths in between 0.9-3.5 G 
are obtained, f or the analyzed e vents. The important conclusion that we extract from 
the analysis by ILin et"aL1 ( 2009 ) is that prominence seismology is possible and works 



well, provided high resolution observations are available. 
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3. Seismology Using Observed Periods and Damping Times 

The damping o f prominence oscillations is a clear feature in many observed events. 
Lin et al.1 ( 2009b show that the amplitudes of the waves passing through two different 



cuts along a thread are notably different. Among the different damping mechanisms 
that have been put forward in order to explain the damping of MHD oscillations in 
prominence plasmas resonant absorption in the Alfven continuum seems a very plausi- 
ble one. The mechanism relies on the non-uniformity of the medium in the transverse 
direction. It was suggested to explain the damping of transverse kink waves in promi- 
nence threads by Arre gui et al.l ( 20081) . The damping of wave modes can be used as an 



additional source of information about the physical properties of promin ence plasmas. 
In the cont ext of transverse l y os cillating coronal loops this was done by lArregui et al. 



(2007) and iGoossens et al] (120081) . Their analytical and numerical inversion schemes 
make use of the simple idea that it is the same magnetic structure, whose equilibrium 
conditions we are interested to assess, that is oscillating with a given period and under- 
going a given damping rate. By computing the kink normal mode frequency and damp- 
ing time as a function of the relevant equilibrium parameters for a one-dimensional 
model, the period, P and damping ratio, P/Td, have the following dependencies 

P P 

P = P(k z , c, l/a, v A i), — = —{k z , c, 1 1 a), (1) 

with VAi the internal Alfven speed, k z the longitudinal wavenumber, c = pi/p e the den- 
sity contrast, and l/a the transverse inhomogeneity length-scale, in units of the tube 
radius, a. In the case of coronal loop oscillations, an estimate for k z can be obtained 
directly from the length of the loop and the fact that the fundamental kink mode wave- 
length is twice this quantity. For filament threads, the wavelength of oscillations needs 
to be measured. Relations Q] indicate that, if no assumption is made on any of the 
physical parameters of interest, we have two observed quantities, period and damp- 
ing time, and three unknowns, density contrast, transverse inhomogeneity length-scale, 
and Alfven speed. There are therefore infinite different equilibrium models that can 
equally well explain the observations. These valid equilibrium models delineate a one- 
dimensional solution curve in the three-dimensional parameter space (c, l/a, VAi)- When 
partially filled threads, with the dense part occupying a length L p shorter than the total 
length of the tube L are cons idered, the period and damping time o f thread oscillations 
are seen to depend on L p /L ( Soler et al. |2010|; lArregui et al.ll201 1 ). Then, one of such 



curves is obtained for each value of the length of the thread. The solutions to the in- 
verse problem are shown in Fig. |2^ for a set of values for L p . It can be appreciated 
that, even if each curve gives an infinite number of solutions, they define a rather con- 
strained range of values for the thread Alfven speed. This figure also shows that the 
ratio Lp/L is a fundamental parameter in order to perform an accurate seismology of 
prominence threads. Because of the insensitiveness of the damping ratio with density 
contrast, for the typically large values of this parameter in prominence plasmas, the ob- 
tained solution curves display an asymptotic behavior for large values of c. This allows 
us to obtain precise estimates for the filament thread Alfven speed and the transverse 
inhomogeneity length scale in that limit for each of the curves. The computation of the 
magnetic field strength from the obtained seismological curve requires the assumption 
of a particular value for either the filament or the coronal density. The resulting curves 
for a typical coronal density and several values of L p /L are shown in Figure |2p. As 



Prominence Seismology 



215 




Figure 2. Determination of prominence Alfven speed (a) and magnetic field 
strength (b) from the computation of periods and damping times for standing kink 
oscillations in two-dimensional prominence thread models and observations of pe- 
riod and damping times in transverse thread oscillations. The observed period and 
damping time are 20 and 60 minutes, respectively, and L — 10 5 km. 



can be seen, precise values of the magnetic field strength cannot be obtained, unless the 
density contrast is accurately known. 



4. Seismology of Flowing and Oscillating Prominence Threads 

The question of whether propagating disturbances detected in chromospheric and coro- 
nal structures are waves or flows has become a hot topic that demands for both obser- 
vational and theoretical assessment in order to perform an accurate seismology. The 
contribution of Hinode to the subject has been of central importance. Flows constitute 
and additional source of information about the physical conditions of oscillating struc- 
tures in conjunction to oscillating periods and damping times. The first seismological 
application of prominence seismology using Hinode observ a tions of flowing and trans- 
versely oscillating threads was presented by Terradas et al.1 (120081) . using observations 
obtained in an active region filament by Oka moto et all (120071) . The observations show 
a number of threads that flow following a path parallel to the photosphere while they 
are oscillating in the vertical direction. The relevance of this particular event is in the 
fact that the coexistence of waves and flows can be firmly established, so there is no 
ambiguity about the wave or flow character of a given d ynamic feature. We have both 
of them in this particular example. lOkamoto et al.l (120071) analyze 6 threads and Table 1 

displays the relevant measured quantities. 

A detailed view of the event analyzed by lOkamoto et all (120071) shows that when a 
given thread is selected and several cuts along its structure are analyzed as a function of 
time, oscillations that are synchronous along the entire length of the thread are found. 
This means that the maxi mum and minimum am plitudes occur at nearly the same time 
for all locations. This led iTerradas et al.l ( 20081) to interpret these oscillations in terms 



of the kink mode of a magnetic flux tube. 

Let us first neglect the presence of mass flows. By using previous theoretical 
results from a normal mode analysis in a two-dimensional piecewise filament thread 
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Table 1 . Summary of geometric and wave p roperties of vertically oscillating flow- 
ing threads analyzed by Oka moto et al.l (l2007h . 2W is the thread length, vo its hor- 
izontal flow velocity, P the oscillatory period, V the oscillatory velocity amplitude, 
and H the height above the photosphere. 



Thread 2W (km) voCkms" 1 ) P (s) ^kms" 1 ) H (km) 



1 


3600 


39 


174 ± 25 


16 


18 300 


2 


16000 


15 


240 ± 30 


15 


12400 


3 


6700 


39 


230 ± 87 


12 


14700 


4 


2200 


46 


180 ± 137 


8 


19 000 


5 


3500 


45 


135 ±21 


9 


14 300 


6 


1700 


25 


250 ± 17 


22 


17 200 



model by iDiaz et all ([2002) and iDvmova & Rudermanl (120051) . iTerradas et all (|2008) 



find that, although it is not possible to univocally determine the physical parameters of 
interest, a one-to-one relation between the thread Alfven speed and the coronal Alfven 
speed can be established. This relation comes in the form of a number of curves relating 
the two Alfven speeds for different values of the length of the magnetic flux tube and the 
density contrast between the filament and coronal plasma. Figure [3] shows the derived 
values by changing the length of the tube from bottom to top and the density contrast, 
from left to right. An interesting property of the obtained curves is that they display 
an asymptotic behavior for large values of the density contrast, typical of filament to 
coronal plasmas, and hence a lower limit for the thread Alfven speed can be obtained. 
Take for instance thread #6. Considering a length of the total magnetic flux tube of 
L = 100 Mm, an overall value between 120 and 350 km s" 1 for the thread Alfven speed 
is obtained. 

Next mass flows are considered. First a simple approximation is made. Consider 
the form in which mass flow along the cylinder, vo, enters in the linear MHD waves 
equations through the differential operator 

d d 

The terms coming fr om the equilibrium flow can, in a first approximation, be ignored 
because, as noted by Id ymova & Ruderma n (2005), inside the cylinder the terms with 



derivatives along the tube are much smaller than those with radial or azimuthal deriva- 
tives. By following this approach the problem reduces to solving a time-dependent 
problem with a varying density profile, p(z, t), representing a dense part moving along 
the tube with the flow speed. After solv ing the two-dimensio nal wave equations one 
finds that the flow velocities measured by Oka moto et all (120071 1 result in slightly shorter 
kink mode periods than the ones derived in the absence of flow. Differences are however 
small and we find period shifts in between 3 and 5%. 

Finally, a more complete approach to the problem has been followed by Terradas et al 



(2008) who consider the numerical solution of the full MHD wave equations, with no 



further approximations. The thin tube approximation is not used, the flow is maintained 
in the equations. Also the density is allowed to change in the simulation as a nonlin- 
ear code is used. The full numerical result confirms the previous approximate results 
regarding the effect of the flow on the obtained periods, and therefore, on the derived 
Alfven speed values. 



Prominence Seismology 



217 



Thread #1 



Thread §2 




800 
_^ BOO 
j= 400 - 

200 - 



1000 2000 3000 
v Ac (km s~') 

Thread #3 




1000 2000 3000 4000 
v Ac (km s~') 

Thread #5 





1000 2000 3000 4000 
v Ac (km s~') 

Thread §4 




1000 2000 3000 4000 
v Ac (km s~') 

Thread §6 




1000 2000 3000 4000 
v Ac (km s~ ) 



1000 2000 3000 4000 
"v Ac (km s~ f ) 



Figure 3. Dependence of the Alfven velocity in the thread as a function o f the 
coronal Alfven velocity for the six threads observed by Oka moto et all (I2007I) . In 
each panel, from bottom to top, the curves correspond to a length of magnetic field 
lines of 100000, 150000, 200000, and 250000 km, respectively. Asterisks, dia- 
monds, triangles, and squares correspond to density ratios of the thread to the coronal 
gas p p /p c ^5,50,100,200. 



5. Summary 

This paper presents recent developments of the MHD seismology technique applied to 
prominence plasmas. The combination of refined theoretical models for prominence 
thread oscillations together with observations allows for the estimation of difficult to 
measure quantities in these objects. Three complementary techniques for the inversion 
of physical parameters have been discussed. They make use of the MHD kink wave 
interpretation for transverse oscillations in filament and active region threads together 
with observational estimates for quantities such as periods, damping rates, and mass 
flow speeds. In general, the solution to the inverse problem is unable to provide a single 
value for all the parameters of interest, since the quantity of unknowns outnumbers that 
of measured wave properties. However, estimates of physical parameters in a narrow 
range of values can be derived. A proper further development of prominence seismol- 
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ogy requires improvements in both theory and observations. Our example using Hinode 
observations demonstrates that high quality observations and proper theoretical analy- 
sis allow flows and waves to become two useful characteristics for our understanding 
of the nature of solar prominences. 
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